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Abstract -Branching random flights are key to describing the evolution of many physical and 
biological systems, ranging from neutron multiplication to gene mutations. When their paths 
evolve in bounded regions, we establish a relation between the properties of trajectories starting 
on the boundary and those starting inside the domain. Within this context, we show that the 
total length travelled by the walker and the number of performed collisions in bounded volumes 
can be assessed by resorting to the Feynman-Kac formalism. Other physical observables related 
to the branching trajectories, such as the survival and escape probability, are derived as well. 



Introduction. — Consider a single walker initially 
emitted from a point source. Once emitted, the walker 
undergoes a sequence of displacements, separated by col- 
lisions with the surrounding medium. When the scat- 
tering centers encountered by the travelling particle are 
spatially uniform, the inter-collision lengths are exponen- 
tially distributed PQ, with mean free path A. The quan- 
tity a = 1/A, proportional to the probability of particle- 
medium interaction along a straight line, is the total cross 
section: its value typically depends on the particle posi- 
tion and speed [T]. At each collision, the incident particle 
disappears, and k particles (the descendants) are emitted 
with probability pk, whose velocities are randomly redis- 
tributed in angle and intensity according to a given prob- 
ability density, which in principle can vary as a function 
of the number k of descendants pQ . Each descendant will 
then behave as the mother particle, and undergo a new 
sequence of displacements and collisions, giving thus rise 
to a branched structure, as illustrated in Fig. [T] Branch- 
ing random walks as described above lie at the heart of 
physical and biological modeling, and are key to the de- 
scription of neutron transport in multiplying media and 
nucleon cascades, spread of epidemics, diffusion of repro- 
ducing bacteria, and mutation-propagation of genes, just 
to name a few ErlH- 



A central question for random walks is to determine the 
occupation statistics of the stochastic paths in a given re- 
gion V of the phase space: for branching or non-branching 
Brownian motion, for instance, this is intimately related 
to such issues as residence times and first-passage proper- 




Fig. 1: A schematic representation of branching random walks 
in a bounded domain, starting on the surface E or inside the 
volume V. For trajectories starting on the surface, the angle 9 
between the normal to the surface and the incoming trajectory 
obeys a density function depending on dimensionality. 



ties [5"Mllj. For exponential flights, the occupation statis- 
tics is naturally defined in terms of two observables: the 
number of occurred collisions in the volume V and the 
total length travelled in F [I]. In principle, one would 
be interested in assessing the full distribution of these two 
quantities: this unfortunately turns out to be a formidable 
task, even for very simple geometries, and analytical re- 
sults are seldom available [T21U5) . 

A somewhat simpler approach consists in deriving for- 
mulas for the moments of the distributions, which are 
sometimes amenable to exact results. For exponential 
flights, important findings have been reported in recent 
years concerning trajectories starting either on the bound- 
ary or inside a given bounded volume, in the special case 
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Pi = 1, i.e., the well-known Pearson random walk 
In particular, it has been shown that the average travelled 
length for walkers starting from the boundary depends 
only on the ratio of the volume over the surface, which is 
a rather counter-intuitive result. In this Letter we extend 
these findings in two directions, based on the Feynman- 
Kac approach [5T]: first, we include more generally ab- 
sorption and branching mechanisms; second, we show that 
other physical observables associated to the process, such 
as escape and survival probabilities, can be easily derived 
within the same formalism. 

Travelled lengths. — Consider a bounded domain 
of non-zero volume V (with finite diameter) and non-zero 
measurable surface E = dV. An exponential flight started 
in V or on its boundary will end when all descendants 
have been either absorbed in V or escaped from E. In or- 
der to keep notation simple, yet retaining the key features 
of the process, we will assume in the following that par- 
ticles evolve at constant speed v, descendants are emitted 
isotropically and independently from each other, and a is 
constant over the volume V. To start with, we denote by 
L m (ro, ujq, t) the m-th moment of the total length trav- 
elled in V by the particle and all its descendants when 
observed up to time t, starting from a single walker in ro 
with direction ujq at time t = 0. It is possible to character- 
ize L m by resorting to the Feynman-Kac formalism [15] . 
which yields the recursive formula 

--L m (r ,uj ,t)^aJ2^rn. J [(L%(r ,t)] 
v at j=a 

+£*L m (r a ,u; , t) + mlyL" 1 " 1 ^, w 0> *), (1) 



for m > 1, starting with L m (ro, Wo, 0) 
L°(ro, u>q, t) = 1 from normalization. Here 



= and 



is the backward transport operator [I], v\ = (k) = ^2 k kpk 
being the average number of secondary particles per colli- 
sion, and fid = 27r d / 2 /r((i/2) the surface of the unit sphere 
in dimension d. Furthermore, we denote by 



(h) n (r Q ,t) = J ^-h(r ,(j} ,t)- 



(3) 



the average over directions, and by ly the marker func- 
tion of the volume V, i.e., ly = 1 when ro belongs 
to V, and ly = elsewhere. Finally, B m j [zi] = 
B m ,j [zi, Z2, ■ ■ ■ iZm—j+i] stand for the Bell's polynomi- 
als [35], Vj = (k(k — l)...(fc — j + 1)) being the falling fac- 
torial moments of the descendant number, starting with 
^o = l. When trajectories are observed up to a time 
t much larger than the characteristic time scale of the 
system dynamics, we can define the stationary moments 
L m (r ,u>o) = lim Z7"(r , u>o, t), provided that the limit 

t—t+oo 

exists. Intuitively, this condition is satisfied when the 



particle losses due to absorptions and leakages from the 
boundaries are larger than the gain due to population 
growth, which is always the case if i^i < 1 [2]. When 
v\ > 1, this typically amounts to further requiring that V 
is below some critical size V c [T|. In the following we will 
assume that V < V c , unless differently specified: the time 
derivative in Eq. Q] then vanishes at large times, and we 
get 

C*L m (r Q , w ) + ml v L m - l (r , w„) 



i=2 



[{L%(v )] =0, 



(4) 



when 



with the boundary conditions L m (ro € E,Wo) 
ujq is directed outward. 

Following [19] . our aim is now to average Eq. 0]over the 
starting position and direction of the walker. As the do- 
main V is bounded, we can safely define the probability 
measures for trajectories born in the domain and for those 
starting on the surface. Choosing the starting coordinates 
uniformly distributed inside V imposes the uniform vol- 
ume probability measure 

dft dV 

where dft is the solid angle element 0- Similarly, an 
isotropic incident flux uniformly distributed on the fron- 
tier E imposes the surface probability measure 
dfl dE 
% ^E 

where ijd = y/n(d— l)r((d— l)/2)/r(d/2) is a dimension- 
dependent normalization constant, equal to twice the in- 
verse of the average height of the d-dimcnsional unit shell, 
and n is the normal entering the surface 0. By means 
of such probability measures, for any function /i(ro,u>o) 
we define its volume average (h}v and its surface average 
(ft.)s, respectively, by 



-(O-n), 



(6) 



(h) v = 
(hh = 



f dr f du> 



due 



(wo-n)/i(r ,«o). (7) 



dE 

~z J n d 

We integrate then Eq. 2] uniformly over all possible initial 
positions and directions (taking into account the isotropy 
property) , and apply the Gauss divergence theorem. This 
yields the recursive formula 



V r 



v )v + a{vi-l)(L m )v 



J=2 



(8) 



1 A /j,— or cquivalcntly (/—randomness in the language of stochas- 
tic geometry I16II23I . 

2 A fA — randomness 1 1 6 1 1231 . The term cos = (f2 ■ n) implies 
that in polar coordinates trajectories starting on the surface must 
enter the domain with density 6 = arcsin(2£ — 1) in two dimensions 
and 6 = l/2arccos(l — 2£) in three dimensions, £ being uniformly 
distributed in (0, 1] ; see Fig. \T\ 
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Equation [5] relates the m-th moment of trajectories start- 
ing on the surface to the different moments (up to order 
m) of trajectories born inside the volume, and as such 
extends to branching random flights the formulas 



V 

■q d — and (IS' 



(L r ' 



m(L)i 



for m > 1 (9) 



previously obtained for Pearson random walks [17H19) . In 
particular, from Eq. [8] the average length (m = 1) reads 



(L) s = Vd -[l + *(v 1 -l)(L) v ], (10) 

which generalizes the celebrated Cauchy's formula (I) s = 
rjd K (also known as the mean chord length property) , orig- 
inally established for random straight lines drawn from the 
surface of the volume [23] and recently shown to rather 
surprisingly apply also to Pearson random flights [F7H20] . 

The term O = (L)s/A is a measure of the opacity of the 
volume, in that it expresses the ratio between the average 
length travelled in V when an isotropic particle flux is im- 
posed at the surface S and the mean free path. Another 
quantity of interest is x = (^^/(VdV/T,), which is the 
ratio between the average length travelled in the actual 
medium V and the length that the particle would elapse if 
V were empty and the paths were straight lines (the mean- 
ing of the denominator stems from Cauchy's formula). 

In general, Eq. [10] depends on the fine details of the 
process and of the geometry, since the term (L)v is not 
universal. However, when the underlying branching pro- 
cess has V\ = 1, Eq. ITUl yields precisely the Cauchy's for- 
mula. In this case, thus, the quantity (L)s would depend 
only on the geometrical ratio V/T, and not on the specific 
details of the random walk. In particular, would be 
independent of the characteristic jump size A. This simple 
property unfortunately does not carry over to higher mo- 
ments of the travelled length. Indeed, for m = 2 we have 



B 



2,2 21,22 



and Eq. [8] then gives 



V 



(L 2 )z = [2(L) v +a (i^-l) (L 2 ) v + ap 2 ((L)U*o))v] ■ 

Bell's polynomials in Eq.[S]are the signature of branching, 
and for m > 2 introduce some extra non- vanishing terms 
(v2 > 0) with respect to Eq. [9]even when v\ = 1. 

In the absence of branching, i.e., when random flights 
can be either scattered or absorbed, with po + pi = 1 and 
Pk = for k > 2, explicit relations for the probability den- 
sity functions of the travelled length can be also derived. 
Under these hypotheses, Eq. [8] reduces to 



(U 



%X \m(IS' 



ao(L r 



(11) 



where oo = Pqc corresponds to the absorption cross sec- 
tion. In presence of absorption, then, (L m )s < VdV/T,, as 
expected. We denote by f(l) and g(r) the probability den- 
sity of the total travelled length for a trajectory started 




Fig. 2: Branching random flights born in the volume V (or 
on the surface E) and traversing a sub-domain V' of V . The 
portion of the travelled trajectories spent in V' is displayed as 
blue dashed lines. 



on the surface or inside the volume, respectively Then, 
Eq. [TTJ can be identically rewritten as 



l m f(l)dl = Vd 



V 



POO 


'TO 




(7 


Jo 


L r J 



r m g(r)dr. (12) 



Integrating Eq. [T2] by parts and using normalization 
J Q g(r)dr = 1 yields the relation 



g(r) = 



V 



f(l)e aol dl 



(13) 



between the two densities f(l) and g(r). In the limit of 
purely diffusive processes (po — > 0), Eq. [T3] reduces to 



g( r ) 



f(l)dl, 



(14) 



a relation originally established for straight paths [24] and 
later extended to Pearson walks [18] . 

Excursions in sub-domains. — Considerable efforts 
have been devoted to the study of the occupation statis- 
tics of some sub-domain V' included in V. This issue has 
been thoroughly investigated, e.g., in the context of resi- 
dence times for Brownian motion (with or without branch- 
ing; see, e.g., [7HTT]) and Pearson walks [19]. Consider a 
branching exponential flight emitted in V: the particle 
and its descendants may enter V', spend some time in- 
side, branch, possibly die in V or escape, then re-enter 
V , and so on, as illustrated in Fig. [2] The total length 
travelled in V' can be straightforwardly assessed by re- 
sorting to the Feynman-Kac formalism mentioned above. 
Indeed, its moments L' m (ro,u;o) satisfy Eq. U] the marker 
function being restricted to the sub-domain V . Then, 
averaging over all angles and positions inside V yields 



Y;[my(L>™- V+ ff (^-l)nv 

(15) 



j=2 
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For m = 1 we have in particular 



By making use of the binomial formula, we finally get 



(L 1 ) 



Vd- 



V 



II 

V 



+ a{v 1 -l)(L'), 



V 



(16) (N m h=Vd^[pi((N + lD v -(N m ) v +po}. (20) 



As a consequence, for trajectories starting on the surface, 
the ratio between the length travelled in V' and that trav- 
elled in V reads 



II 
V 



1 +jrQ>i -l)(L') v 
l + afa- 1)(L) V 



(17) 



This ratio generally depends on the geometry as well as 
on the walk features. However, for branching flights with 
v\ = 1, we have (L')s/(L)s = V'/V, i.e., we recover 
the elegant ergodic-type property that applies to Pearson 
walks [inj. 

Number of collision events. — The stationary m-th 
moment N m of the number of collisions in V performed by 
a branching exponential flight starting from ro in direction 
ujq can be also assessed by resorting to the Feynman-Kac 
formalism: as shown in [15] . N m satisfies 

m 

£*iV m (r , w ) + aj^^mj [(N%] + 



i=2 



i—k 



fc=l ^ ' 3=0 



0. 



(18) 



which closely resembles Eq.|U At the boundaries, we have 
N m (ro <G £,a>o) =0 when u>q is directed outward. Then, 
by averaging Eq. [18] over starting positions and directions 
in V we get 



{N r 



Vd- 



V 



m 
fe=l 



m—k 

E 



Vj{B m - k4 [(JV*)n(r )]>v 



+ - 1) (N m ) v + y E,"i(B m , i [(N%(v )])v 

j=2 



which generalizes the results previously found for Pearson 
walks |25| . For the average collision number we have 

(N)v = v d j;a[l + (v 1 -l)(N)v]. (19) 

Now, from exponential flights being a Markovian process it 
follows that (N)v = <?{L)v [H|- Hence, we have also O = 
cr(L)s = (N)s, which amounts to saying that the opacity 
can be expressed in terms of the mean number of collisions 
in V. Similarly as done for the lengths, the number of 
collisions in a sub-domain V' can again be computed by 
using lv in Eq. [T5] 

For non-branching flights, Vj = for j > 2 and, from 
[zi\ = Zj and Bj t o [z,] = for j > 1, we have 



V 
Vd^a 



Pi 



E 

k=Q 



(N 1 



i — k\ 



(N r] 



Po 



In the absence of branching, it is also possible to explic- 
itly derive the collision probabilities. We denote by f(i) 
and g(j) the collision number probability for a trajectory 
started on the surface or inside the volume, respectively. 
Then, Eq. [50] can be identically rewritten as 



i=0 



V 
Vd^a 



Ebi(i 

3=0 



] m ]gU)+Po 



By equating the terms of the series, and imposing that the 
relation holds for arbitrary m > 1, we get then 

f(j) = nd\o [ Pl9 (j - i) - gU)+PoS it i] , (21) 

Sij being the Kronecker delta. Resumming over j yields 
in particular /(0) = 1 — r}d^crg(0). For Pearson walks, 
Po — > and we obtain 



V 

fti) = Vd^cr [g(j - 1) - g(j)} 



(22) 



Escape probability. — When the volume V is 
bounded, the probability i?(ro, u>o) that a trajectory start- 
ing from ro in direction u>o never visits the exterior of V 
satisfies [15] 

- o; • V ro i?(r , u; n ) + aR(r , w ) = cxl v G [{R) n \ , (23) 

where G[z] = J2kP kzk ^ s * ne generating function as- 
sociated to the descendant number distribution pk and 
i?(ro £ S, Wo) = when u>o is directed outward. By 
definition, the quantity Q — 1 — R represents the escape 
probability from the volume. Then, averaging Eq. !23l over 
all initial positions and directions yields 

(Rh = r,d^[(G[{R}n])v-(R)v]. (24) 

The terms (R)v and and have a simple probabilis- 

tic meaning, namely, the probability that a particle born 
uniformly and isotropically in the volume V , or entering 
the body isotropically from the boundary, respectively, 
is absorbed with all its descendants in V [T]. Similarly, 
(i?)n(r ) represents the probability that a particle born 
isotropically at ro is absorbed (with all its descendants) in 
V. Equation [M] relates therefore the spatial behavior of 
the probability R to the probabilities pk, via G[z\. 

When walkers can not be absorbed in the domain (po = 
0), trajectories must necessarily escape from the bound- 
aries, and we have (R)s = 0. This rather intuitive result 
can be understood as follows: developing Eq. 1241 and using 
the normalization ^2^ Pk = 1 yields 

V °° 

(Rh = %^E^ [« R )a)v - (R)v] ■ (25) 

k=2 
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Since (-R)q < (R)n {(R)n is a probability), we immedi- 
ately get that (-R)s < 0: then, the probability that a 
particle entering the body isotropically from the bound- 
ary is absorbed with all its descendants in V must vanish, 
as expected. The same applies to (R)v 

In the absence of branching, G\z\ = p a + p\z with p + 
pi = 1, and we get 



(Rh = Vd^ao [i-{R)v], 



a d-dimensional generalization of a theorem originally de- 
rived for purely absorbing media [T] and extended to dif- 
fusive and absorbing media in three dimensions in |18j . 

Survival probability. — A fundamental quantity 
for random walks is the survival probability <St(ro,u;o), 
namely, the probability that at time t at least one particle 
is still in V. When V < V c , for long times ^(ro, u; ) — > 0, 
which physically means that the combined effects of ab- 
sorptions in V and leakages through E are sufficient to 
compensate the population growth due to branching (if 
any), and trajectories almost surely go to extinction. How- 
ever, when v\ > 1 and V > V c , branching paths have a 
finite probability of surviving indefinitely in V , and there 
exists a non-trivial limit St(ro,u>o) S(ro,u>o) > 0, 
depending on the starting position and direction. For 
branching Brownian motion, for instance, finding the 
asymptotic survival probability is a long standing is- 
sue [25][37]- For branching exponential flights in bounded 
domains, when V > V c the probability of ultimate survival 
S(r ,uJo) satisfies [15] 

-<a - W ro S(T ,oj )+aS(r ,u: ) = aF[(S}u], (27) 

where F[z] = ET=i a kZ k , with a k = (-l) fe+1 ^/fc!. At 
the boundaries, S must vanish when u>o is directed to- 
wards the exterior of V. Averaging Eq. [27] over all initial 
positions and directions yields then 

(S)v = r,d^<T[{F[(S) a ]) v -{S)v]. (28) 

Unfortunately, the complex nature of the alternating series 
in F[z] seems to prevent from drawing general conclusions 
based on Eq. [2E] 

Perspectives. — The approach proposed in this 
Letter based on the Feynman-Kac formalism allows 
the volume-averaged properties of branching exponential 
flights to be derived by relying upon a minimal number 
of simplifying hypotheses. The results presented here are 
fairly general, and as such apply to a broad class of physi- 
cal and biological systems. Moreover, no assumptions are 
made concerning the shape of the domain, and our findings 
remain valid for non-convex bodies, including for instance 
domains with holes. Further investigations are ongoing in 
three directions: including the effects of anisotropies and 
heterogeneities, both in the source terms and at the colli- 
sion events [28], allowing for other (possibly unbounded) 



jump distributions, such as Levy Flights |29| . and taking 
into account reflective or mixed boundary conditions [19] . 



The authors wish to thank S. Mohamed for carefully 
checking formulas by Monte Carlo simulation. 
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